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THE ATTITUDE DETERMINATION OF A THREE 
AXES STABILIZED SPACECRAFT 

by 

Ai Chun Fang 

Telemetry Computation Branch 
Information Processing Division 
Tracking and Data Systems Directorate 

ABSTRACT 

This report presents an approach for attitude determination of a three axis sta- 
bilized spacecraft. It describes the employment of two independent vectors whose com- 
ponents with respect to two reference coordinate axes can be known. Such two coordinate 
axes are right-handed three mutually orthogonal axes systems of which one is fixed in 
the spacecraft and the other has known orientation in space. The relationship between 
these two systems are then expressed into mathematical equations with the use of those 
two reference vectors. By computing the derived equations, the attitude of the space- 
craft can therefore be obtained. 
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THE ATTITUDE DETERMINATION OF A THREE 
AXES STABILIZED SPACECRAFT 
by 

Ai Chun Fang 

Telemetry Computation Branch 
Information Processing Division 
Tracking and Data Systems Directorate 


INTRODUCTION 

Assume that there is a first coordinate system fixed in the spacecraft and the spacecraft is 
stabilized about all three axes. Also assume that a second reference coordinate axes is being 
used whose orientation in space is always known. The attitude control system of the spacecraft 
is required to provide the stabilization and to maintain the alignment of the first system to the 
second coordinate axes. Any deviation of the first system from the second one is detected by 
the sensing instruments on-board the spacecraft. The complete sensor data are then employed 
to determine two independent vectors in space, such as the sun's position vector and the earth s 
magnetic field. These vectors are applied to obtain the attitude of the spacecraft. This report 
presents the analysis of the attitude determination of such a spacecraft and the formulation of 
required mathematical equations. The application of this approach is also discussed and shown 
in simple example. 


C oordinate Axes Definitions 

For convenience, it is necessary to make the following definitions: 

(a) The Spacecraft Body Axes.-The spacecraft body coordinate axes x, y, z is a set of right 
handed mutually orthogonal axes fixed in the spacecraft with their origin at the center of mass. It 
is referred as first coordinate system mentioned above. 

(b) The desired Orientation System.-This system X, Y, Z with its known orientation in 
space is called the second reference coordinate axes sometimes. The orientation of this orthogonal 
axes system is always the desired one for the spacecraft body axes. That is, the x, y, z axes are 
to align with X, Y, Z respectively if there is no attitude error. 
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(c) Orbital Coordinate Axes. -The orbital coordinate axes X 0 , Y 0 , Z 0 is used only for the 
spacecraft with circular orbit. This system has its center at the instantaneous location o 0 of the 
spacecraft with 

X 0 = spacecraft's velocity vector 

Y 0 = normal to the orbital plane 

Z g = the downward local vertical defined as the unit vector from the location of the spacecraft 
to the center of the earth. 

(d) The inertial coordinate system is an earth center right handed orthogonal system such 
that the Z. axis coincides with the north pole and the X i Y plane is the equatorial plane in which 
X 4 axis passes through the vernal equinox. 

The Attitude Determination Equations 

It is assumed that the direction cosines of the reference vectors in the desired system X, Y, 

Z are known and that in the x, y, z axes have been determined. The attitude of the spacecraft can 
thus be obtained by utilizing the information of these reference vectors. 

Suppose that the transformation taking the desired orientation system into the body axes is 

x - A n X + a 12 Y + A, 3 Z 

? = A 21 x + A 22 y + A 23 z (1) 

z = a 31 X + a 32 Y + a 33 z 

Let v be a reference vector in space and let 



be the components of V in the body axes and in the desired coordinate system respectively. Then 

V-x = Ajj V-X + A 12 V-Y + A 13 V-Z 

V-y = A 21 v-x + a 22 v-y + a 23 v-z 
v-* = a 31 v-x + a 32 v-y + a 33 v-z 
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This means 


V „= A U V 1 + A ,2 V 2+ A 13 V 3 

V y = A 21 V 1 + A 22 V 2 + A 23 V 3 ^ ^ 

V , = A 31 V 1 + A 32 V 2 + A 33 V 3 

Equation (3) implies that any vector v in the X, Y, Z system written in column form, may be trans 
formed into x, y, z system by the application of the matrix [A] . . That is 


V 


~v,l 

X 


1 

V 

y 

< 

11 

V 2 

V 

z 


V 3 

L_ 


where 


(4) 


A 

A in 

A, , 

li 

1 2 

13 

A 

A 

A . 

n 21 

22 

23 

A 

A 

A., 

31 

32 

33 


(5) 


The matrix [A] is defined by the coefficient A. . (i, j = 1, 2, 3) in equation (1) and is now to be 
determined. 


Since [A] is a 3 x 3 matrix, it requires three reference vectors for obtaining the solution of 
A. . Let 5, E, and F be such required three vectors. Suppose that [D] d , [E] d and [F] d are 
computed and also that [D] b , [El b and [F] b have been obtained. By applying equation (4) to the 
vectors D, E, and F respectively, it yields 


and 


D 


r D ri 

X 


i 

D 

= [A] 

D 

y 

2 

D 


D 3 

Z 


L 3 J 


E 


E i~] 

X 


1 

E 

y 

= [A] 

E 2 

E 

Z 


A 


1 

A 


i — 

A 

F 

y 

< 

II 

F 2 

A 

L 


L F 3_ 


(6) 


(7) 


( 8 ) 
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Combining these three column matrix equations into a single square matrix equation, gives 


~D 

X 

E 

X 

F ~ 

X 


" D i 

E ! 

p r 

D 

y 

E 

y 

F 

y 

= [A] 

° 2 

£ 2 

f 2 

D 

Z 

E 

t 

F 


_°3 

E 3 

•*) 

1 


It may be written in this simple notation 


where 


[B] = [A] [Q] 


[B] = 


F 

X 

F 


and 


[Q] = 


F 2 ' 


°3 E 3 


( 9 ) 


(10) 


If the determinant of [Q] is not zero, then [Q] is non- singular and possesses an inverse [Q] -1 . 
By multiplying the matrix equation (10) with [Q]' 1 , gives 

[B] [Q ]" 1 = [A] [Q] [Q ]- 1 


or 


[A] = [B] [Q]* 1 (11) 

It should be noticed that knowing two reference vectors D and E the third vector F can be 
defined from them. For simplicity, let D and E be unit vectors. By defining 


F = Dx E 


the column matrix [FL and [FI are 

a b 


[F] 


” F 1 


D 2 E 3 - D 3 E 2 

F, 

_ 

D, E, - D, E, 

2 


3 1 13 

A 


5 e 2 - d 2 e l 


( 12 ) 
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(13) 


F n 


D 

E 

- D 

E 

X 


y 

Z 

Z 

y 

F 

_ 

D 

E 

- D 

E 

y 


Z 

X 

X 

Z 

F 


D 

E 

- D 

E 

L *j 


X 

y 

y 

X 


[FL 


If 77 is the angle between vectors D and E, the angle 77 is calculated from 

cos 7 ? = Ej + D 2 E 2 + D 3 E 3 

Then the matrix [Q] can be written as 


E . D 2 E 3- D 3 E 2 


[Q] = ! D 2 E 2 D 3 Ej - B 1 E 3 


D 3 E 3 D , E 2- D 2 E 1 


and its determinant is 


If 7 j / 0, then 


det [Q] - | Q | = sin 2 


iQl ^ 0 

The matrix [Q] is non- singular and has an inverse [Q ]' 1 where 


(14) 


(15) 


[QT 1 


sin i) 


D 3 - Ej cos 77 D 2 - E 2 cos 77 D 3 - E 3 cos 7/ 
Ej - D x cos 77 E 2 - D 2 cos 77 E 3 - D 3 cos 77 


D 2 E 3 - °3 E 2 °3 E ! - D ! E 3 D . E 2 " °2 E ! 


(16) 


Now equation (11) becomes 


~ A u 

A 12 

r 

<f 


X 

E 

X 

D 

y 

E 

Z 

- D E ^ 
x y 


Dj - Ej cos ?? 

D 2 -E 2 cos 77 

D 3 

- E 3 

COS 77 

A 21 

A 22 

A 23 

1 

s in 2 77 

D 

y 

E 

y 

D 

Z 

E 

X 

- D E 

X Z 


Ej - D 3 cos 77 

E 2 - D 2 cos 7) 

E 3 

- D 3 

COS 7] 

A 31 

A 32 

A 33 


D 
L 1 

E 

Z 

D 

X 

E 

y 

- D E 
y x 


D 2 E 3 - D 3 E 2 

D 3 E 1 - D ! E 3 

D 1 

E 2- 

D 2 E 1 


(17) 


This directly shows 


A n = {(D Dj + E x Ej) - (D E t + ^ E % ) cos r, + (D y E, - D E y ) (D a E s .- D 3 E 2 )}/sin 2 7, 

A 12 = °2 + E X E 2> - < D x E 2 + c 2 E x > cos V + (D y E - D E y ) (» 3 E, - ^ E 3 )}/sin 2 7, 

A i3 = ((°x °3 + E X E s) - < D X E 3 + c 3 E x) cos V + (D y E - D E y ) (Dj E 2 - D 2 E^l/Sin 2 77 
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22 


*32 


33 


= <(D y D 5 + E y E t ) - (D y + D t E y ) cos r, + (D t E x - D x E,) (D 2 E 3 - D 3 E 2 )}/sin 2 
= { (D y D 2 + E y E 2 ) - (D y E 2 + D 2 E y ) cos v + (D z E x - D x E,) (D 3 E J . D J E 3 )}/sin 2 
= {(D y D 3 + Ey E 3 ) - (D y E 3 + D 3 Ey) cos V + (D, E x - D x E ) (Dj E 2 - D 2 E^J/sin 2 

= { (° 2 D j + E x Ej) - (D, Ej + \\ E ) cos Tj + (D x E y -D y E x ) (D 2 E 3 _ D 3 E 2 )}/sin 2 

= { ( D , D 2 + E * E j> ~ (D E 2 + D 2 E 2 ) cos V + (D x E y - D y E x ) (D 3 E ; - ^ E 3 )}/sin 2 

= { {D °3 + E E a) - <D E 3 + D 3 E ) cos v + (D X Ey - D y E x ) (D 3 E 2 - D 2 E^/sin 2 


•n 

V 

V 

V 

V 

V 


(18) 


From the above equations the deviation of the body axes from the desired orientation system can 
therefore be known. For FORTRAN program, equation (17) is recommended. 


Euler Angle Trans formation 

The attitude error of the body axes x, y, z from the desired orientation system X, Y, Z can 
be represented by a sequence of three Euler angle rotations. The first is that a rotation about 
axis p through angle <p, the second about axis q through angle , and the third about axis r through 
angle 8 . Thus there are twelve different sequence of Euler angle rotations that can be formed such 
as 121, 123, 131, etc. However, any one of them is available for choosing to represent the space- 
craft attitude. For instance, let the sequence 312 be selected. That is, yaw angle 0 around the 
Z axis, roll angle 4> about the displaced X axis, and pitch angle 8 about the doubly displaced Y 
axis (see figure 1). The notation to be used for the matrix which performs each of these three 
Euler rotation is 


R (0) = 


COS l p 
- s in i p 


0 


sin i|/ 0 

cos i// 0 

0 1 


1 


R x «>) = 


o 

0 


0 0 

cos <fi s in <^> 
- s in 4> cos <j> 


v*> 


cos 8 0 -sin 6 

0 10 
s in Q 0 cos 8 
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Therefore 


V 

X 


V 

V 

y 

= R y (0) \( 4 >) R,« 

V 2 

V 

z 


_ V 3 


( 19 ) 


By comparing equation (19) with equation (4) it gives 


or 


[A] = R y (9) R x (<£) R z (4>) 



cos 8 

0 

- sin 6 


1 

0 

» ' 


cos p 

sin 

0 

[A] 3 

0 

1 

0 


0 

cos <p 

sin d> 


- sin i/i 

cos \p 

0 


s in 8 

0 

cos 6 


0 

L_ 

- s in 4> 

cos <p 


0 

0 

1 


cos cos i/i - sin 8 sin<£sinyi 
- cos 4> s in \p 

sinfi cosi/i-cos 9 sin ^ sin ^ 


cos Q sin i p + sin 9 s in <j!> cos p 
cos <p COS l p 

sinfi’sinV'-cos 9 sin cp cos 4 1 


- s in 9 cos <p 
sin p 

cos 6 cos <p 


( 20 ) 


The roll angle 4>, pitch angle 9, and yaw angle p can thus be derived 


and 


Hence 


A 13 = - sin 9 cos 4> 
A 23 = sin <p 


A = cos 6 cos <p 


A 21 = - sin f cos <ii 


<P - sin -1 (A 23 ) 


6 - tan" 


of the angle of 4> is limited to - 90 < 4> < 90 then cos <P is always positive and 


sin p - - 


c ° s * r v^tal 


( 21 ) 

( 22 ) 
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Figure 1. Attitude Error Representation: Roll 4 >, Pitch 8 , Yaw 


Thus 

<P = sin 1 (- 

Application 


There are three reference vectors whose components in the desired orientation system can 
always be known. They are the sun position vector, the earth magnetic field, and the downward 
local vertical. In order that Equation (17) or (18) can be applied, the direction cosines of each 
reference vector in the body axes must also be found. At present, three kinds of sensors have 
been developed for this purpose. The sun sensor measuring the direction cosines of the sun posi- 
tion vector, the earth sensor providing the information of the downward local vertical, and the 
magnetometer giving the components of earth magnetic field. Since most spacecraft contains only 
two kinds of sensor because of limited space, the given information can thus provide two available 
reference vectors to be used instead of three. 

It should be noticed that the desired orientation system is different from one spacecraft to 
the other. The transformation for calculating the components of the reference vector in that system 
must be formulated according to the requirement of each spacecraft. 
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For example if a spacecraft is in a circular orbit with an on-board magnetometer, the 
components of the earth magnetic field M in the body axes will be measured. Let M x , M y , and M x 
be such three components. Then 


and 


3 = m/|m| 



"d ~ 

X 


i 

[D] b = 

D 

y 


M 

y 




M 

X 


(24) 


The spacecraft position unit vector R is known from orbit determination and will be indicated 
in orbit 3A tape. If R ( , R 2 , and R 3 are the three direction cosines of vector R in the inertial co- 
ordinate system at the time when the magnetometer makes its measurements, then 


cos a = Rj //Rf + R 2 


s i n a = R 2 //r* + R 2 
cos v = R 3 

sin v = /r2 + r2 


(25) 


Since the instantaneous position O q of the spacecraft is known, a local coordinate system X e , Y e , 
Z e with its origin at O 0 can be defined where 

X e is pointing toward north 
Y e is pointing toward east 

Ze = Z fl 


At the point O 0 the direction cosines (M N , M E , My) of the rpagnetic field M in the X e , Y e , Z e axes 
can be derived from the actual knowledge of the field. Suppose now the inertial coordinate system 
is the desired orientation for the spacecraft body axes. The (M, , M 2 , M 3 ) of M in that system should 
be calculated by using the following coordinate transformation 




- cos V 

COS O' 

- s in a 

- s in v cos a 


X 

m 2 

= 

- COS V 

sin cr 

COS O' 

- sin v sin cr 


m e 



sin 

V 

0 

- COS V 

J 




(26) 


9 



Hence 



( 27 ) 


Also assume that the spacecraft has on-board earth sensor providing the information for com- 
puting the components L x , L y , and L t of the downward local vertical Z 0 . Then 


and 



Knowing ID] b , [D] d , [E] b ,and [E] d , equation (18) can be computed immediately. 

Next suppose that desired orientation of the spacecraft is changed from the inertial coordinate 
system to the orbital coordinate axes because of the purpose for special experiment. The column 
matrices [D] b , [E] b will remain in the same form, but [D] d and [E) d are different. In this 
case 



where i is the orbit inclination (i.e. angle between orbit plane and equatorial plane), and 
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Figure 2. The X^, Y £ and Z,£ Axes in the Local Coordinate System 



Substituting these values into equation (17) or (18) the problem of attitude determination is then 
solved. 


CONCLUSIONS 

The approach described here may well be suited for attitude determination of non- spinning 
spacecraft consisting of any two of the three on-board instruments; namely, the sun sensor, the 
magnetometer, and the earth sensor for providing the informations of two reference vectors. The 
only restriction for this approach is that the two selected reference vectors must not be parallel. 

The idea of this approach was first practiced on Delta- PAC in Jan. 1969 when the author was 
assigned to formulate the attitude determination equations and error prediction formulae for Delta- 
PAC project. The derived equation were exercised with actual flight data obtained from Delta- PAC 
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after it was launched on Aug. 9, 1969. The results determined the attitude behavior of Delta- PAC 
and evaluated its control system. However, the Delta- PAC attitude determination equations were 
limited to be used only for that particular spacecraft and cannot be directly applied to others if 
not modified. Therefore, the more generalized attitude determination equations for a three axes 
stabilized spacecraft seem to be needed and are contained in this report which are expected to be 
applicable for most non- spinning spacecraft. 
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